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Abstract 

We have studied the problem of all sky search in reference to continuous gravita- 
tional wave particularly for such sources whose wave-form are known in advance. We 
have made an analysis of the number of templates required for matched filter analysis 
as applicable to these sources. We have employed the concept of fitting factor (FF); 
treating the source location as the parameters of the signal manifold and have studied 
the matching of the signal with templates corresponding to different source locations. 
We have investigated the variation of FF with source location and have noticed a sym- 
metry in template parameters, 6t and (pT- It has been found that the two different 
template values in source location, each in 0t and 4>t, have same FF. We have also 
computed the number of templates required assuming the noise power spectral density 
Sn{f) to be flat. It is observed that higher FF requires exponentially increasing large 
number of templates. 



1 Introduction 

Gravitational wave (GW) Laser Interferometer antennas are essentially omni - directional 
with their response better than 50% of the average over 75% of the whole sky (Grishchuk 
et al., 2000). Hence the data analysis systems will have to carry out all sky searches for 
its sources. We know that the amplitude of intense GW believed bathing the earth is very 
small, as compared to the sensitivity of GW detectors, and is further masked by the dom- 
inant noise. In these circumstances, continuous gravitational wave (CGW) sources are of 

* e-mail : dcsrivastava @ now-india. com 
^e-mail: ssahay@iucaa.ernet.in 

^Present address: Inter University Centre for Astronomy and Astrophysics, Post Bag 4, Ganeshkhind, 
Pune-41 1007, India 



1 



prime importance because for such sources we can achieve enhanced signal-to-noise ra- 
tio (SNR) by investigating longer observation data set. However, a long observation time 
introduces modulation effects, arising due to the relative motion of the detector and the 
source. As a consequence, there results redistribution of power in the forest of side bands 
resulting into the reduction of the expected power due to amplitude modulation (AM). The 
problem of all sky search gains another dimension in view of the fact that there are reasons 
to believe the presence of intense GW sources whose locations and even frequencies are 
not known. Amongst such sources pulsars occupy an important position. Similar to all sky 
search one will also have to do all frequency search. All sky all frequency search is the holy 
grail of gravitation pulsar astronomy. In this paper we confine ourselves to the problem of 
all sky search. 

Search of CGW without a priori knowledge appears to be computationally quite de- 
manding even by the standard computers expected to be available in the near future. For 
example, in the case of bandwidth 10^ Hz, observation time 10*^ sec. and star's minimum 
decay time of 100 years one would require 10^^ T flops computer (Frasca, 2000). Very fast 
computer and large memories with ample amount of disk space seems inevitable. How- 
ever, a choice of optimal data processing and a clever programming is also integral part 
of a solution to this problem. Amongst these the pre-correction of time series due to the 
Doppler modulation before the data is processed may be a method, which will reduce the 
computational requirements. In reference to this, Schutz (1991) has introduced the concept 
of patch in the sky as the region of space throughout which the required Doppler correction 
remains the same. He has also demonstrated that the number of patches required for 10^ 
sec. observation data set and one KHz signal would be about 1.3 x 10^^. However, the size 
of the patch would also depend on the data analysis technique being employed. 

Matched filtering is the most suitable technique for the detection of signals from sources 
viz., pulsars whose wave form is known. The wave forms are used to construct a bank of 
templates, which represent the expected signal wave form with all possible ranges of its 
parameters. The time of arrival, source location, frequency of the signal, ellipticity of the 
source and its spin down represent important parameters of GW emitted by a pulsar. For 
detection of GW we check if the cross correlation of the templates with the corresponding 
data set exceeds the preassigned threshold. We introduce in the next section the criterion 
of the fitting factor (FF) (Apostolatos, 1995) applicable to such analysis. We consider 
the source location as parameters of the signal manifold and investigate the matching of 
the waveforms corresponding to different source locations. In section 3 we compute the 
number of templates required for all sky search. A discussion of the results is provided in 
the section 4. 

2 Matched filter analysis: Templates 

The bank of templates is matched, in practice, to only a discrete set of signals from among 
the continuum of possible signals. Consequently, it is natural that all the signals will not 
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get detected with equal efficiency. However, it is possible to choose judiciously the set of 
templates so that all the signals of a given amplitude are detected with a given minimum 
detection loss. The standard measure for deciding what class of wave form is good enough 
is the FF. It quantitatively describes the closeness of the true signals to the template mani- 
fold in terms of the reduction of SNR arising due to the cross correlation of a signal outside 
the manifold with the best matching templates lying inside the manifold. If the FF of a 
template family is unity the signal lies in the manifold. If the FF is less than unity the 
signal lies outside manifold. 

Even if the signal discrete templates lie within the template manifold it would be un- 
likely that any of the actual templates used would correspond to the signal. The parameters 
describing the search template (source location, ellipticity, etc.) can vary continuously 
through out a finite range of values. The set of templates characterised by the continuously 
varying parameters is of-course infinite. However, in practice the interferometer output 
must be cross correlated with a finite subset of the templates whose parameter values vary 
in discrete steps from one template to the next. This subset ("the discrete template family") 
has measure zero on the manifold of the full set of possible templates ("the continuous 
template family"), so the template which most closely match a signal will generally lie in 
between the signal and the nearest of the discrete template family. The mismatch between 
the signal and the nearest of the discrete templates will cause some reduction in SNR. This 
would mean that the members of the discrete template family must be chosen so as to render 
acceptable loss of SNR. 

The study of templates has been made by many research workers in time domain 
[Schutz (1991), Krolak (1997), Brady et. al. (1998), Brady and Creighton (2000), Jara- 
nowski et al. (1998) and Jaranowski and Krolak (2000)]. However, the analysis in the 
frequency domain has the advantage of incorporating interferometer's spectral noise den- 
sity. In order to determine the number of templates required for matched filtering analysis, 
we make use of the formula for FF as given by Apostolatos (1995) i.e. 



where h{f) and /ir(/; ^t, 4>t) represent respectively the FTs of the actual signal wave form 
and the templates. 

In earlier papers we have made an analysis of continuous gravitational wave output 
through Fourier transform (FT) (Srivastava and Sahay, 2002 a, b). These papers are referred 
in the sequel respectively as I and II. It has been established there that the AM of CGW data 
output results into redistribution of power at four additional frequencies / ± 2frot, f =t frot 
in accordance with the frequency modulation (FM). Hence it is sufficient for the analysis of 
FF to consider only the frequency modulated FT. The results obtained in paper- II regarding 
FT of frequency modulated data output i.e. (11-27, 28) may be arranged, making use of the 
symmetry property of the Bessel functions, as follows. 



FF 



max 
0,<l> 



{h(J)\hAf:OT-C>T)) 



(1) 




Hf) 



[{sin(7^ + Q) - sin(7e + Q - vQ} + 
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i {cos(7^ + Q) - cos(7^ + Q - uQ}] + 

JoiZ) E ^^[(3^W-A'V)-*(A'W + 3^V)] + 



m=l 
fc=oo m=cxD 



(2) 



fc=l m=— oo 



A* = sin(7^+ Q-m7r/2) 

3^ = cos(7e +Q- rmi/2) 

U = sin uC^o cos m(^o ~ ^) ~ v {'^'-'^ "^('Co — 5) — sin m6} 

V = cos z/^o cos m(^o — 6) + — sin z/^o sin m(^o — 5) — cos m5 




Now it is straight-forward to compute FF. To understand the procedure let us consider a 
source at location {9, (p) = (25", 30") emitting frequency fo = 0.5 Hz. In order to evaluate 
the summation given in Eq. @ let us note that the value of Bessel function decreases 
rapidly as its order exceed the argument. Accordingly, for the present case it is sufficient to 
take the ranges of k and m as 1 to 1610 and —4 to 4 respectively. We wish to analyze the 
data set for To = one sidereal day. To maximize FF over 9 and 0, we first maximize over 
(j) by fixing (pr to some arbitrarily selected value say, (pT = <P = 30°. Having done this we 
maximize over 9 by varying 9t in discrete steps over its entire range i.e. 0" to 180°. The 
results obtained are plotted in Fig. (p. It is remarked that in order to compute the inner 
product of two waveform hi and /12 which is defined as 



where * denotes complex conjugation,~denotes the Fourier transform of the quantity under- 
neath (a(/) = a{t)exp{—2nift)dt) and Sn{f) is the spectral density of the detector's 
noise. One would require to integrate the expression over the band width of Doppler mod- 
ulated signal. The band width may be determined by computing the maximum value of 
the Doppler shift. In accordance with Eq. (1-23) the Doppler shift is ~ 10^^ fo- However, 
for the present case we consider the band width equal to 0.002 Hz. In a similar manner 
one may fix the 9t and obtain the variation of FF with template parameter (px- Figures (g) 
and (|3|) represent such a plot for a signal of fo = 25 Hz emitted by a source located at 
9 = r&nd(p = 35° and 220°, respectively. 

The following points in reference to these plots may be noted. 



(i) The FF is unity for 9t = 25°, 155° [Fig for 0t = 35°, 145° [Fig @] and for 
(Pt = 220°, 320° [Fig d)]. 





(4) 
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(ii) It is found that the dependence of FF on the template variables 9t and 0t may be 
expressed via 



pp = g-0.00788(e - Ot? (5) 
PP ^ g-0.01778(</. - </.t)2 



(iii) The oscillatory behaviour is more or less typical in waveform that match well or bad 
depending on their parameters. However, we are content with the technique we have 
employed as the region of such artificial facets fall in the region of the FF < 0.25. 

Finally, we conclude this section by noting the symmetry property of the FF with template 
parameters. A closer look of the graphs and the remark (i) above reveal the following 
symmetry property. The FF is symmetrical under following transformations. 

Ot — >7r-eT 0<0t <7T (7) 

4>T ^ TT — < < TI" (8) 

(pT — > Stt — (f)T TT < 0r < 27r (9) 

Let us note that these symmetry properties are based on our results for one day observation 
time. The generic nature of the symmetries may be established only after studying the 
variation of FF with Tg. Unfortunately, we could not make this analysis because of our 
limitations on the computational facilities. 



3 Number of Templates 

It is important to study the problem of number of templates for all sky search in the light of 
FF. The results of the previous section reveal that the grid spacing A9 in the ^—parameter 
of templates may be expressed symbolically as a function of FF, fo,To,9 and (p as 

Ae = j^{FF,ue,<p,T,) (10) 

Similarly, we have 

Ac^ = g{FF,ue,(p,n) (11) 

In the literature it is reported that the choice of grid spacing Acf) depends insignificantly 
on (f), whereas A9 depends as sin26' (Brady and Creighton, 2000). However, in order to 
arrive at such conclusions one may have to make analysis for different values of 9 and 0. 
Unfortunately, due to the limited memory and the efficiency of the computer we could not 
make a study of this aspect. 



5 



In view of this, Eqs. (Q) and (|^) may be equivalently expressed as 



J^(FF,0.5,25^30°,T) 



(0.00788)"Mn(FF) 



1/2 



(12) 



g{FF,25,r,35°,T) 



(0.01778)"Mn(FF) 



1/2 



(13) 



For any chosen value of FF one can determine AO and A0. However, there is no unique 
choice for it. Our interest would be in the assignment of A9 and Acj) such that the spacing 
is maximum resulting into the least number of templates. As we have mentioned earlier, 
there is stringent requirement on reducing computer time. Accordingly, there is serious 
need of adopting some procedure/formalism to achieve this. For example, one may adopt 
the method of hierarchical search given by Mohanty and Dhurandhar (1996) and Mohanty 
(1998). This search is carried out in two steps. At the first level one would start with 
template bank with a coarse spacing in the parameter space but with a lower threshold. 
At the next level a more finely spaced set of templates and a higher threshold would be 
used but only around those templates of the previous level which exceeded the previous 
threshold. 

However, an important issue related to the problem of number of templates is regarding 
the study of the behaviour of number of templates with FF for different fo and Tq. We 
have made an investigation of this aspect. We assume a source location (9, cp) = (1°, 30°). 
We choose some value of FF, say 0.995. Now we maximize FF over by taking 
4> = 4>T = 30" and choose the spacing AO so as to yield the selected FF. In the case 
under investigation AO is found to equal 4.5 x 10^^. Thereafter we maximize over by in- 
troducing spacing Acj) in the so obtained bank of templates and determine the resulting FF. 
The results obtained may be expressed in the form of a graph such as shown in Figs. (Q) 
and (|^). Interestingly, the nature of these curves are similar. We have obtained a best fit to 
the graphs as 

A^TempZates = exp[a - 6x + cx^ - (ia;^ + ex''] ; 0.85 < x < 0.99 (14) 

where a, 6, c, d and e are constants as given in Table (|1]). 

Let us note from the graphs, for sake of comparison, that the number of templates 
required for FF equal to 0.97 are respectively 1.44 x 10^°, 3.5 x 10^° and 5.5 x 10^° 
for observation data sets of 30, 120 and 365 days and fo = 50 Hz. Similarly to analyse 
the observation data set of 120 days of GW frequencies 20, 50 and 100 Hz the number of 
templates required are 1.22 x 10^°, 2.16 x 10^° and 5 x 10^°. It is observed that higher FF 
requires exponentially increasing large number of templates. 

4 Discussion 

In view of the complexity of the FT, which contains trignometric as well as Bessel func- 
tions; one has to be careful in computing FF. We have found useful to employ the 
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fo 


T 


a 


b 

xlO""^ 


c 

X 1 


d 

xin~^ 


e 

xin~^ 




30 


2138.05 


2071.43 


7225.73 


6239.43 


2036.14 


50 


180 


2317.05 


-71.3155 


1746.61 


2146.55 


944.931 




365 


2382.96 


216.917 


2464.42 


2464.42 


1031.54 


20 




2047.55 


-794.473 


3564.56 


4650.68 


1945.87 


50 


120 


2266.59 


4269.44 


15655.0 


17509.5 


6484.51 


100 




2360.23 


-206.906 


1158.27 


1491.01 


733.520 



Table 1 : Coefficients of the best fit graphs obtained for the number of templates 
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Romberg integration using Pade approximation. We have used (i) QROMO of numeri- 
cal recipes instead of QROMB as the former takes care of singularities, and (ii) RATINT 
routine for Pade approximation. 

We have noticed marked symmetries in all sky search in both 9 and space for one day 
observation time. It has been found that the two different template values in source location, 
each in 9t and 0t, have same FF. Accordingly, computation burden will be reduced by a 
factor of four. However, it is not clear whether this symmetry property can be established 
analytically as well. The source location, because of these symmetries is uncertain and 
some other analysis is to be adopted for getting the exact location. We have computed the 
number of templates assuming the noise power spectral density Snif) to be flat which is 
justified as the bandwidth is extremely narrow. 

The issues of optimum template parameterization and placement, and the related com- 
putational burden have been discussed in the literature by several authors notably by Sathyaprakash 
and Dhurandhar (1991), Dhurandhar and Sathyaprakash (1994), Owen (1996), Apostolatos 
(1995, 1996), Mohanty and Dhurandhar (1996), Mohanty (1998), Owen and Sathayaprakash 
(1999). The question of possible efficient interpolated representation of the correlators is a 
problem of current interest and remains still unsolved. 
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Figure 2: Variation of FF with 0^- 
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Figure 4: Variation of number of templates with FF for fixed at different To. 
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Figure 5: Variation of number of templates with FF for fixed To and of frequencies. 
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